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ABSTRACT 

We investigate the impact of spatial variations in the ionized fraction during reionization on tem¬ 
perature anisotropies in the CMB. We combine simulations of large scale structure to describe the 
underlying density field with an analytic model based on extended Press-Schechter theory to track 
the reionization process. We find that the power spectrum of the induced CMB anisotropies depends 
sensitively on the character of the reionization epoch. Models that differ in the extent of the “patchy 
phase” could be distinguished by future experiments such as the Atacama Cosmology Telescope (ACT) 
and the South Pole Telescope (SPT). In our models, the patchy signal peaks at I ~ 2000, where it 
can be four times larger than the kinetic Sunyaev-Zel’dovich (kSZ)/Ostriker-Vishniac (OV) signal 
(ATtot — 2.&^K). On scales beyond I ~ 4000 the total Doppler signal is dominated by kSZ/OV, but 
the patchy signal can contribute up to 30% to the power spectrum. The effect of patchy reioniza¬ 
tion is largest on scales where the primordial CMB anisotropies dominate. Ignoring this contribution 
could lead to significant biases in the determination of cosmological parameters derived from CMB 
temperature measurements. Improvements in the theoretical modeling of the reionization epoch will 
become increasingly important to interpret the results of upcoming experiments. 

Subject headings: cosmology: theory - cosmic microwave background - large-scale structure 


1. INTRODUCTION 

Cosmic Microwave Background (CMB) anisotropy ex¬ 
periments have now constraine d temperature fluctu¬ 
ations down to sca les of 4' llR.eadhead et al.l 12004 
iHolzanfel et aUHoOPll resulting in a much improved un¬ 
derstanding of inhomogeneities during the time of de¬ 
coupling at z ~ 1100 and of the global properties of the 
Universe. In combination with data from Supernovae la, 
measurements of the local expansion rate, galaxy cluster¬ 
ing studies, and observations of the Lyman a forest, the 
CMB data has lead to the establishment of a “standard” 
cosmological model. We live in a Universe dominated 
by dark energy and with significantly more dark matter 
than baryons. Structure grew from a scale independent 
spectrum of primordial Gaussian fluctuations, in good 
agreement with predictions of inflationary models. 

In the CMB community, theoretical and experimental 
interest is shifting to the study of secondary anisotropies 
on smaller angular scales. These are caused by fluctua¬ 
tions in the distribution of baryons and dark matter in 
the redshift regime z ~ 0 — 30. Future data on the dif¬ 
ferent secondary effects will constrain the way structure 
formation proceeded from the linear into the non-linear 
regimes. Secondary anisotropies can be divided into 
three categories: gravitational leasing effects, inverse 
Compton scattering of CMB photons by hot plasma, and 
Doppler related anisotropies. 

Mass concentrations along the line of sight such as 
clusters, sheets, and filaments lead to weak deflection 
of the CMB photons, or Gravitational Lensinq (for a 
review see chapter 9 of iBartelmann fc Schneid'CTl2001fl . 
Gravitational leasing of the GMB can be used as a cos- 
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mological probe in various ways. The leasing effect on 
the CMB power spectrum could he lp break degeneracies 
between cosmological param eters llMetcalf fc SiHdlTQQSl: 
iStomoor fc EfstathioiJ ITo^ . Leasing induces depar¬ 
tures from Gaussianity on GMB maps that can be used 
to reconstruct the spatial di s tribution of the leasing mass 
llZaldarriaga fc Seliaklll994 iHii fc Okamotoll2002ll . The 
GMB leasing signal can be correlated with galaxy leas¬ 
ing shear providing additional information. Thus, leas¬ 
ing can be used to probe the evolution of gravitational 
clustering, constrain the properties of the dark energy, 
the shape of the matter power spectrum, and neutrino 
physics (see e.g. IVan Waerbeke et all 1200(1: IHiJ I2002t 
iKaplinghat et al 

Inverse Compton scattering in the hot intracluster 
medi um, also called the thermal Sunyaev-Zel’dovich ef¬ 
fect (iZebdovich fc SimvaevlIlQO^ . changes the spectrum 
of the CMB photons, leading to cold spots (decrements) 
in the microwave background at frequencies below 217 
GHz and to hot spots (increments) at frequencies above. 
The effect is proportional to the line-of-sight integrated 
pressure. Because it is independent of redshift, it is a 
unique probe of collapsing structures out to z ~ 3 (see 
Figure [ 7 ||. The thermal SZ effect has been measured 
in foll ow-up observations of X-ray clusters (se e for in¬ 
stance iBirkinsha.wlll 9991: iCarlstrnm et al .1 l2flflnfl . Gom- 
parison of the cluster SZ temperatures with the X-ray 
measured temperatures leads to constraints on the an¬ 
gular diameter distance. The SZ effect also leaves a sig¬ 
nature in the CMB power spectrum on small scales where 
the primordial CMB vanishes. BIMA (the Berkeley Illi¬ 
nois Maryland Association) and CBI (the Cosmic Back¬ 
ground Imager) cl aim detections of this effect at > 2a 
signif icance levels l|Readhead et al.ll2004 IHolzanfel et alJ 
[1993). CBI infers a LSS clustering amplitude on 8 Mpc/h 
scales of as = 0.9, which is at the upper level of con¬ 
straints from cluster observations and weak lensing (for 
a recent compilation of experimental results see table 5 of 
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iTegmark et al.ll2ti(ill . Future large angular scale SZ sur¬ 
veys such as SPT promise to measure the cluster abun¬ 
dance as a function of mass and redshift, which will offer 
the possibility of constraining the matter density and the 
equation of state of the dark energy. 

Finally, Doppler related effects are produced by the 
scattering of CMB photons off electrons moving as a re¬ 
sult of the structure formation process. This process, also 
known as the kinetic Sunya ev-Zerdovich effect when ap - 
plied to clusters of galaxies llSunvaev fc ZerdovichllQSill . 
leads to hot or cold spots, depending on whether the ion¬ 
ized baryons move toward or away from the observer. 
The frequency dependence of the photons is left un¬ 
changed, except for tiny relativistic corrections. These 
“Doppler” induced anisotropies are the only known way 
to measure the high redshift large scale velocity field. 

Although our simulations account for each of these ef¬ 
fects, our focus in this paper will be on the part of the 
Doppler effect generated during the epoch of reionization. 
This epoch is currently not well understood observation- 
ally. 

The first year release of the WMAP data hinted at a 
large optical depth owin g to reionization, Tri = 0.17 ± 
0.05 llSnergel et alJ 1200,111. through mea surement of the 
“reionization bumD’ MZaldarriagal llQQ TT in the temper - 
ature polarization cross correlation lIKogut et aUl^OO.lll . 
WMAP places only a weak constraint on the relative 
height of the acoustic oscillation peaks, so these data 
alone are not sufficient for constraining Tri, because it 
is strongly degenerate with the tilt of scalar fluctuations 
ris and the baryon density parameter uj}, = In 

combination with smaller scale CMB measurements, the 
data can be described by homogeneous reionization of 
the universe taking place at redshift z ~ 14 ± 3, or a 
Tho mson scattering optical depth Tri — 0.13 ± 0.02 (see 
e.g. iReadhead et al.ll2004 [Bond et al.ll2003ll. This esti¬ 
mate is also robust to the addition of other datasets to 
the analysis, such as the SDSS galaxy clustering survey 
llTegmark et ahll^OOdH {Tri = 0-124lQ;g57 ) or the Lyman 
alpha forest llSeliak et al.ll2004ll {Tri = 0-133to;o45)- 

Since the large scale polarization measurements are 
sensitive only to the line-of-sight integrated free electron 
density, they do not require the time evolution of the ion¬ 
ized fraction to be a simple step function. Observations 
of the evolution of the Lya optical depth in the absorp¬ 
tion spectra of the highest redshift quasars at z ~ 6 
when combined with the CMB constraints favor a more 
complicated ionization history. The observations show a 
distinct Gunn-Peterson trough l|Cunn fc Petersonlll96,5(l 
and point to a rapidly evolving neut ral fraction indica¬ 
tive of the final stages of reionizati on llBecker et an i2001t 


iFan et al.ll2004jl . Re¬ 
cently, there has also been an interpretation of the rela¬ 
tively high temperature of the Lya forest at z ~ 2 — 4 
as evidence of an order unity change in the ionized frac - 
tion at z < 10 llTheuns et al.l200'^lHui fc Haimanl2003jl . 
althoug h this depends on th e properties of He II reion¬ 
ization llSokasian et alJl2002ll . 

Numerical simulations of the star formation epoch have 
shown that a first generation of low metallicity (Popula¬ 
tion III) sources can survive negative feedback from the 
UV background th ey produce if shielded in massive ha - 
los for some time l|Scliaeredl2003t iBromm et al.ll2002jl . 


If they are heavy, these stars could produce an order of 
magnitude more ionizing photons per baryons than “nor¬ 
mal” stars. Besides a more complex temporal behavior, 
the ionization fraction could have spatial variations; it 
could be “patchy” for a period of time while reionization 
is not complete. Naturally, we would like to have an ob¬ 
servable that sheds light on the duration and structure 
of the reionization epoch. 

The shape of the large scale polarization signal in 
the CMB could be used to constrain the evolution of 
the global ionization fraction, without resorting to infor¬ 
mation from secon dary anisotropies l|ZaldarriagallI997t 
IHu fc Holder! 1^. However, these scales are affected 
by sample variance limitations, so the available informa¬ 
tion is rather limited. Moreover, before such an approach 
can be succ essful, issues of foreground removal need to 
be cla rified llBennett et al.ll2?)0^lde Oliveira-Costa et alJ 
l2?inl . In any event, large scale polarization measure¬ 
ments will not yield information about the morphology 
of reionization, as this is confined to sub-degree scales. 

Several analytical models for patchy reionization have 
been presented in the literature. They vary significantly 
in the assumed size and time evolution of spherically 
shaped bu bbles, and whether the ionized regi ons are 
correlated jKuioxgl^lJ 11994 iSantos et alJl200,^ or not 
(fCruzinov fc HiJlT99^ Some models employ a prescrip¬ 
tion for correlating the ionizing sources by using the bias 
that can be calculated for the dark matter halos in which 

All authors 


they presumably reside 
agree that on scales below 4', where the primordial CMB 
signal falls rapidly owing to photon diffusion, the Doppler 
effect induced by patchiness could contribute enough to 
use it as a tool to study the reionization epoch. 

Modeling patchy reionization more accurately requires 
the use of numerical simulations to incorporate effects of 
non-linear clustering and radiative transfer. The mor¬ 
phology of the ionized regions depends sensitively on the 
star formatio n model (see e.g. the reionization study 
companion bv iBarkana fc Lo^l200Ill . source properties, 
feedback processes, and radiative transfer (these effects 
have for example been modeled inlAbel fc Wandeldl20r)4 


Gnedin fc ,Taffell200H iCiardi et al.11200,31: iSokasian et alJ 


2mi is. These simulations have given us evi¬ 
dence that reionization was much less homog e neous 
than previously thou ght llMiralda-Escude et all II998I: 
iBarkana fc LoebI 1200111 . Reionizat ion tends to proceed 
from high to low density regions l|Sokasian et alJ 1^00,31 
l2?iol iCia.rdi et al.l l200.3|l and recombinations seem to 
play a role subdominant to large scale bias. Stromgren 
spheres of neighboring protogalaxies overlap, and over¬ 
densities will harbor large ionized regions. Because HII 
regions extend to larger radii than the correlation length 
of galaxies, building an analytic model solely on local 
galaxy properties appears to be difficult. Besides un¬ 
certainties over the physics in current radiative transfer 
calculations of reionization, memory and CPU require¬ 
ments pose a serious limitation, so that these simulations 
have so far only been performed as a post-processing step 
on scales of up to 1 0 Mpc/h, corresponding to ~ 6 ar- 
cminutes on the sky llFurlanetto et al.ll2004a^ . If the HII 
regions are of comparable size, sample variance becomes 
a problem for these simulations. 

In this paper, we use a hybrid approach between an 
analytic model for the formation of HII regions, and 
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smoothed particle hydrodynamics (SPH) simulations of 
large scale structure. The advantage over analytic mod¬ 
els of patchy reionization is that our sources follow the 
complex clustering behavior of dark matter and baryons. 
An advantage over currently feasible full radiative trans¬ 
fer calculations of the reionization epoch is that we can 
investigate the morphological properties of reionization 
on scales an order of magnitude larger with rather small 
memory and CPU requirements. Our prescription does 
not directly address many of the (uncertain) physical 
details of reionization related to star formation, feed¬ 
back processes, dumpiness, recombinations, and radia¬ 
tive transfer. However, by combining these properties 
into a single parameter, we are able to explore the ba¬ 
sic parameter dependencies, and we believe it is a good 
starting point to help us understand the morphological 
properties of cosmological reionization and its influence 
on the CMB. 

We will describe the analytic model for the morphol¬ 
ogy of p atchy reionization suggested bv lFiirlanetto et id] 
(l2004btl in the next section. In this description, the size 
distribution of HII regions is derived in a way analogous 
to the Press-Schechter formalism for the halo mass func¬ 
tion, so that sources are biased towards large scale matter 
concentrations. The model has only one free parameter, 
the efficiency of the first sources to ionize the surrounding 
medium. 

In Section |3 we describe our implementation of the 
model into SPH simulations as a post-processing step. 
Since the costs in memory and CPU consumption in this 
scheme are relatively low, one can imagine testing a wide 
array of behaviors for the ionization efficiency of the hrst 
sources against experiments. 

One of the parameterizations we employ, in Section 
0 attempts to do justice to the apparently contradict¬ 
ing measurements of a large integrated Thompson opti¬ 
cal depth in CMB polarization measurements that indi¬ 
cates an onset of reionization at redshifts beyond z = 
15, and a rapid decline of the neutral fraction around 
z = 6 — 7, seen in measurements of the Gunn-Peterson 
optical depth by measuring spectra of distant quasars. 
We model this by assuming two succeeding reionization 
epochs, in which HII regions first are produced by low 
metallicity (Pop. HI) stars, and later expand further 
during an epoch dominated by Pop. H stars, so that on 
the whole patchiness lasts longer. 

In SectionElwe investigate whether the next generation 
of ground based bolometric arrays will be able to measure 
reionization and distinguish between different scenarios. 
The speci fications of the A tacama Cosmology Telescope 
(ACT) |22113) and the South Pole Telescope 

(SPT)^ llRuhl et al. 1200411 will be used to predict how 
well these experiments can measure CMB temperature 
power spectra in the region where secondary effects dom¬ 
inate over primordial CMB anisotropies. We hnd that 
ACT/SPT should be able to distinguish between sudden 
homogeneous reionization and patchy reionization at a 
high level of significance, even if the total optical depth 
in both scenarios were the same. 

Since the patchy signal peaks on scales where the pri¬ 
mordial anisotropy dominates, cosmological parameter 

see http://www.hep.upenn.edu~angelica/act/act.html 

® see http://astro.uchicago.edu/spt/ 


estimation may become biased. We will show in Section 
15. 21 th at precision experiments such as Planck should take 
this bias into account either by focusing on their polar¬ 
ization data or by adding a reionization parameter to 
their analysis. 

Our numerical simulations of large scale structure and 
the power spectra we used to generate realizations of the 
primordial CMB assume a cosmology in agreement with 
WMAP constraints: ildm = 0.26, Ha = 0.7, Hf, = 0.04, 
fluctuations on 8Mpc/h scales of erg = 0.9, and no tilt 
(ns = 1) or running (as = 0) of the spectral index. 


2. ANALYTIC MODEL FOR PATCHY REIONIZATION 

We review the model for the growth of HII reg i ons dur - 
ing reionization proposed bv iFurlanetto et ffiT ll2004bll . 
The standard way to describe ionization of the IGM is 
to associate an HII region with a single galaxy. The size 
distribution of ionized regions then follows directly from 
the halo mass function when the ansatz niion = (,‘^gai is 
made, where nigai is the mass in a collapsed object. The 
parameter C is the efficiency factor for ionization, for ex¬ 
ample composed as C = where fesc is the 

escape fraction of ionizing photons from the object, /* is 
the star formation efficiency, the number of ioniz¬ 

ing photons produced per baryon by stars, and Urec is 
the typical number of times a hydrogen atom has recom¬ 
bined. The efficiency factor is a rough combination of 
uncertain source properties, but encapsulating a variety 
of reionization scenarios in it can be regarded as a start¬ 
ing point to gain insights in the morphological properties 
of the partly ionized phase. 

According to t he extended Press-Schechter model 
(iBond et alJllQQlll . the collapsed fraction (or the frac¬ 
tion of baryons that lie in galaxies) in a region of size r 
depends on the mean overdensity of that region, Sr, as 


/co//(PImin) — erfc 


Sc(z) - Sr 

y^2[a^(rmin) “ ^^(r)] 


( 1 ) 


Here, (r) is the linear theory rms fluctuation on scale r 
and rinin is taken to be the radius that encloses the mass 
niinin (at average density p) corresponding to a virial tem¬ 
perature of lO'^K, at which atomic hydrogen line cooling 
becomes efficient and Sc(z) is the numerical factor 1.686 
scaled to today using linear theo ry. The redshift depen- 
dence of the minimum mass fsee lBarkana, fo I:oeUl2flfll1l 
can be described by 

~ 3.0 X 10®(1 + z)-i'^M0 . (2) 


The mass in galaxies that can create enough luminous 
sources to fully ionize all hydrogen atoms is inversely 
proportional to the ionizing efficiency, so one requires 
that 

fcoll > • (3) 


Hence, we can define a barrier which fluctuations have 
to cross for their baryonic content to become ionized: 

Sr > Sj;(m,z) = dc(2:)-V2erfc“^(C“^)[cr^(r'i„in)-CT^(r)]^/^ . 

( 4 ) 

Because it assumes Gaussian fluctuations on the mass 
scale m, this formalism can be applied only to mass scales 
larger than the typical size of collapsed objects. 
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3. SIMULATION OF PATCHY REIONIZATION 

Our approach is to use large scale simulations of 
the cosmic web as a basis for applying the model 
described in the last section as a post-processing 
step to generate HII regions during the reioniza¬ 
tion phase. The underlying large scale structure 
was simula ted using the parallel Tree-PM/ SPH solver 
GADGET (ISoringel. Yoshida fc Whit^ 1200111 , based on 
the fully conservative im plementation of SPH by 
ISoringel fc HernauistI ll2002ll . Here, we used the results 
of runs with a boxsize of 100Mpc//i and 216^ particles, 
with parameters corr e spond ing to the G-series runs of 
ISoringel fc HernauistI (1200311 . In what follows, we em¬ 
ploy a simulation that included only “adiabatic” gas 
physics; i.e. the gas can heat or cool adiabatically and 
be shocked, but we do not include radiative effects or the 
consequences of star formation and associated feedback 
processes. Snapshots of the simulation are produced ev¬ 
ery light crossing time interval. This leads to 77 outputs 
between z = 0 and z = 20. These s i mulat ions are de¬ 
scribed in more detail in lWhite et alJ (1200211 . 

To generate inhomogeneously ionized regions accord¬ 
ing to the model described in Section [3 we discretize 
the matter fluctuations in our simulation boxes into 256^ 
cells. The overdensity, <5, is smoothed with a top hat win¬ 
dow function. 

To find the smoothing scale at which a given cell is ion¬ 
ized we have to keep in mind that it could also be ionized 
by photons originating from neighboring regions. Thus 
we should smooth the density on all possible scales to see 
whether a given point was above the ionization threshold 
of equation for some smoothing scale. In practice, we 
start at large radii (comparable to the simulation box 
size) and record the smoothed overdensity as we smooth 
logarithmically on progressively smaller scales. When 
a cell first crosses the barrier (which depends only on 
redshift and ionization efficiency), it is deemed ionized. 
If later, as a function of decreasing scale, it crosses the 
barrier downwards, this means that the region had been 
initially ionized by a neighboring overdensity. 

Figure^shows the random walk of dp, the density fluc¬ 
tuation scaled to the present using linear theory, for three 
different regions inside our box (at redshift z = 14.2) 
with increasing (T^(to) (decreasing smoothing radius r). 
It also shows the ionization efficiency dependent barrier, 
from equation ( 0 , as the solid curve. The short-dashed 
curve describes a region of high overdensity that self- 
ionizes. In the dotted curve at a^(rn) ~ 2.2 the barrier 
is crossed downwards, so the volume element was ionized 
by sources in neighboring cells. The long-dashed curve 
corresponds to an element that did not ionize at this 
redshift. 

To sample a wide range of bubble sizes, and obtain 
a smoothly varying ionization fraction, we varied the 
barrier in radial direction throughout each box by ad¬ 
justing the redshift dependent minimal ionization radius. 
The reionization information for each of the 256^ cells is 
stored. A 8Mpc/h deep cut through the stacked outputs 
between z ~ 12 and z ~ 17 is shown in Figure 13 for 
one of the models we will describe in the next section, in 
which the ionization efficiency is constant, C, = 60. The 
periodic boundary conditions of our simulation boxes are 
apparent in this plot. 



Fig. 1.— The behavior of the overdensity scaled to today, 5o, 
for four regions from our 256^ calculation. The abscissa in this 
plot is the rms fluctuation, dependent on the applied smoothing 
scale. The solid curve is the barrier given from equation 13). The 
region corresponding to the short dashed curve crosses the barrier 
for ionization at ~ 1.8. The region described by the dotted 

curve wanders below the barrier again, but has been ionized by 
sources in a neighboring region with higher density. Finally the 
long dashed curve represents a region in which the gas stays neutral 
at this particular redshift {z = 14.2). 


To compute the electron scattering and gravitational 
leasing effects from the content of the simulation boxes, 
we only need to store two-dimensional maps of the prod¬ 
uct of pressure and volume (Compton scattering), the 
velocity weighted free electron density (Doppler effect), 
and the matter density (leasing). We simulate regions 
on the sky with an angular extent of 1° each, at a resolu¬ 
tion of 256^ pixels. In the small angle approximation we 
can simply project the content of each simulation box 
onto a plane in its center. At high redshifts, z > 5.7, 
the field of view exceeds the boxsize, so we use the peri¬ 
odic boundary conditions to cover it. During the patchy 
regime, only those gas particles that are located in an ion¬ 
ized cell according to the information we stored earlier 
contribute to the signal. At later times all particles are 
assumed ionized. We translate and rotate the positions 
and velocities of the particles randomly before doing the 
projection onto the plane. The gas particle properties are 
distributed over surrounding pixels weighted by the SPH 
kernel. For dark matter particle masses we use the cloud- 
in-cell algorithm. Photons are traced through the planes 
from a regular grid close to the observer towards the last 
scattering surface. We produced Gaussian random fields 
of the primordial GMB using spherical harmonic decom¬ 
position coefficien ts generated with the public ly available 
code CMBFAST (ISeliak fc ZaldarriagaiilQQlf 

In models of patchy reionization, GMB anisotropies 
are caused by two types of contributions induced by 
Doppler scatterings: perturbations in the baryon den¬ 
sity Wb, given by 5iOh = Pb/ptot — 1, and local changes 
in the ionization fraction 5xe- These produce a change 

® http://www.cmbfast.org 
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z=16.61 z=15.53 z=14.19 z=13.16 z=12.24 


Fig. 2.— Time development between redshifts z = 12 and 2 : = 17 of an 8Mpc/h deep slice through an inhomogeneously ionized (^ = 60 
(constant) box, with increasing ionization fraction. The outputs with periodic boundary conditions are simply stacked behind each other 
here. In order to obtain realistic maps, the individual outputs where randomly rotated and translated in our simulations. At the low 
redshift end, the bubbles become comparable in size to the box. 


in temperature of the CMB blackbody. The total con¬ 
tribution to the temperature anisotropy is given by the 
integral over conformal time ^ 

^r^(n) = -—npivo) [ d? 7 [a" 2 e"'"”(’'^Xe(? 7 )]n • q 

J-CMB c J 

(5) 

with 

g = (l + ^)(l + <5a;b)v. (6) 

Xe{r]) 

The Thomson scattering optical depth is ut, Tri is the 
optical depth from the observer to conformal time rj, and 
n is the line of sight unit vector. 

To improve our understanding of reionization, various 
parameterizations of C and its derivatives may be com¬ 
pared to data from the next generation of experiments. 
For each redshift, the only free parameter in our model 
of patchy reionization is the efficiency of the first sources 
which ionize the surrounding intergalactic medium. In 
contrast to radiative transfer simulations, our implemen¬ 
tation is not limited by CPU or memory related prob¬ 
lems. For each box we need to Fourier transform ~ 50 
times back and forward to do the smoothing. This takes 
~ 100 minutes per model on a Xeon 3.2 GHz computer. 
The implementation may be parallelized easily, hence the 
realization of a large number of parameterizations of is 
feasible. 

Even in the context of our simplified model, there are 
two additional limitations to the accuracy of our present 
simulations. First, because the mean overdensity on 
scales larger than the simulation box is artificially set to 
zero in numerical simulations, in order to achieve peri¬ 
odic boundary conditions, we ex pect a systematic bias in 
the overall ionization fraction. iBarkana fc LoebI ll2004fl 
estimated this effect and showed that it should be less 
than 1% in simulations of 100 Mpc/h size. A different 
bias arises from our finite mass resolution. Our lack of 
structure on scales smaller than Thoa:/(Vpart.)^^^ leads 
to a slight delay in the onset of reionization. Very small 
overdensities on scales below 0.5 Mpc/h, that could har¬ 
bor ionizing sources with HII regions around them, are 
not captured by our analysis. 

^ During homogeneous reionization and at lower redshifts, where 
galaxy clusters are present, this is traditionally called the kinetic 
Sunyaev-ZePdovich effect, a t higher redshifts the Ostriker-Vishniac 
JQstnker_&_^ishni^^98^ effect 


4. RESULTS FOR DIFFERENT TIME DEPENDENCE OF 
THE IONIZATION EFFICIENCY 

In our analysis we will compare three different reion¬ 
ization scenarios. These are tuned such that they all 
lead to a comparable integrated optical depth, — 
0.125 — 0.133. The different histories of the fraction of 
ionized volume elements Q{z) are shown in FigureEl 

Model A describes a universe that undergoes ho¬ 
mogeneous and instantaneou s reionization at redshift 14 . 
This is the standa rd scenario l|Seliak fc ZaldarriagallQQOl: 
iTewis et al.ll200(Hl assumed in likelihood analyses of large 
scale CMB polarizati on anisotropy (in i ts correlation 
with temperature, see iKogut et alJl200,i: ISoergel et alJ 
l2f)0l . This model has r^i = 0.133. 

Model B is the patchy model described above, with 
the assumption of constant source properties. An ion¬ 
ization efficiency (/ = 60 leads to reionization beginning 
at z = 19 and concluding by z = 12. This model has 
Tri = 0.130. 

Model C exhibits extended patchiness. We describe 
this by assuming a first generation of metal free sources 
starting at redshift 20 with an ionization efficiency of 
C = 200 (it has been suggested that their photon output 
cou ld be 10-20 times higher than that of normal stars, 
see lSchaeredl2003l iBromm et al.ll200^ LSomerville et alJ 
1200311 . In our model, these sources sustain themselves 
only to redshifts around 15, because their hard pho- 
tons dissociate fragile they need for co oling (e.g. 
iMackev et al.l200^lYoshida et al.l2003Ll2004^ . However, 
the first sources are assumed to leave pockets of ionized 
medium which begin to harbor Pop H stars. The initial 
mass function becomes less top heavy with time, and we 
describe the resulting Population H stars with an effi¬ 
ciency ^ = 12 (see, e.g. Furlanetto et al. 2004c). This 
leads to a total optical depth of Tri = 0.125. Our model 
has a monotonically evolving ionization fraction. Double 
reionization scenarios with an intermittent decrease of 
th e ion i zation fra ction have for instance been suggested 
in iCenI ll2003ll and iWvit he fc Lo^ ll2003l ) but seeni diffi¬ 
cult to realize in practice llFurlanetto fc Loebll2004ll . We 
note that another way of reconciling CMB with quasar 
spectra observations is by resorting to recombinations 
that leave 0.1% of the univ erse neutral until z = 6 — 7 
llSokasian et al.ll2003il2003l . 

In FigureEl we plot for all three models the fraction of 
ionized volume elements as a function of redshift, Q{z). 









































6 


0.6 

'n' 

O' 

0.4 


0.2 


0 I-^^^^^.. . I 

6 8 10 12 14 16 18 20 


Fig. 3.— Evolution of the fraction of ionized volume elements 
Q{z) in our calculations of reionization. Solid curve corresponds 
to instantaneous reionization at redshift 2 = 14, yielding an op¬ 
tical depth of Tri ~ 0.125 for standard cosmological parameters. 
The dashed curve represents our patchy model in the case of a con¬ 
stant ionizing efficiency of sources. The dotted curve is for a model 
where a first generation of metal-free stars can survive the negative 
feedback of H 2 photodissociation for a limited period of time, then 
formation of these stars comes to a halt, leaving a network of HII 
regions. After a while. Population II stars are born inside these 
regions and their photons gradually lead to a homogeneously ion¬ 
ized universe (Model C). Model C spends the longest time in the 
partly ionized regime and will therefore lead to the largest signal 
of the three models which have comparable total optical depth. 


In Model A, the neutral fraction drops instantaneously 
at z = 14, in Model B sources of a constant ionization 
efficiency C = 60 lead to a brief partly ionized phase, and 
in Model C the ionization fraction “freezes in” for some 
time, while metal free Pop. Ill stars cease to exist and 
leave residual ionized bubbles for “normal” Pop. II stars 
to be created and to ionize the medium fully at z = 7 — 8. 

We divide the total Doppler signal into three differ¬ 
ent redshift regimes in Figure 2| The maps have a side 
length of 1° and are smoothed with a Gaussian beam of 
width 6 ~ 1', corresponding to a multipole number of 
Imax — 10000, comparable to the angular resolution of 
ACT, and slightly below that of SPT. It is expected that 
during the interval from x = 0 to 2 = 3, more strongly 
clustered regions make the Doppler effect maps highly 
non-Gaussian. We show this by the solid line (PDF) 
in Figure El In the picture in the center (z = 3 — 11) 
structure formation is less advanced. The correspond¬ 
ing histogram can be well-approximated by a Gaussian, 
as shown by the long-dashed curve in Figure El Sig¬ 
nals created by radially moving inhomogeneities during 
this epoch are referred to as t he Ostriker-Vishniac ef¬ 
fect llOstriker fc Vishniadll98^ . The plot on the right 
of Figure 01 is the patchy epoch of reionization Model 
B. Again, this epoch has a largely Gaussian morphology 
(short-dashed line of Figure EJ, because structures on 
scales of the ionized bubbles are just beginning to col¬ 
lapse. According to this figure, large bubbles form on 
scales of tens of Mpc. The redshift regimes represented 
by the left and middle picture have the same properties 


in all three of our reionization scenarios. 

The upper panel of Figure El shows that in the homo¬ 
geneously reionizing universe, the largest portion of the 
Doppler signal comes from low redshifts, the kSZ regime. 
The upper panel represents the case of homogeneous 
reionization. Model A. Gontrary to the thermal SZ ef¬ 
fect, the signal still carries information about higher red- 
shifts, z < 10. We can use this principle to gain knowl¬ 
edge about the details of reionization taking place at high 
redshifts. The additional contributions from larger red- 
shifts (z > 10) in Model G are shown in the lower panel 
of the figure. 

For comparison, we show in Figure[71how different red¬ 
shift regimes contribute to the thermal SZ effect (where 
AT/T = —2y for low frequencies). The plot shows that 
this signal is basically saturated at z = 3; the thermal 
SZ effect receives most of its contribution from galaxy 
clusters in the more nearby universe. 

5. OBSERVABILITY OF PATCHY REIONIZATION WITH 
FUTURE CMB EXPERIMENTS 

In this section, we assess the observability of patchy 
reionization with future experiments. ACT and SPT 
should be able to distinguish different reionization sce¬ 
narios with high significance by measuring temperature 
power spectra. In their cosmological parameter analyses, 
all-sky experiments such as Planck will have to account 
for patchy reionization as a possible source of bias and 
need to rely on their polarization data in order to obtain 
proper constraints. 

5.1. Power spectral constraints from ACT and SPT 

The effect of patchy reionization on the CMB power 
spectrum is of similar magnitude to that of the Doppler 
effect induced by variations in the mean density. The first 
source population is heavily biased so the patchy reion¬ 
ization signal peaks on larger angular scales. Figure El 
shows that for the regime in which the ionization frac¬ 
tion is roughly 50% (which is where most of the signal 
comes from) the bubbles reach sizes of tens of comov¬ 
ing Mpc. In Figure El we plot the different contributions 
to the Doppler power spectrum. Patchy reionization in 
Model C with ^ = 12 — 200 lasts from redshifts 20 to 
7. The signal imprinted in the CMB from this epoch is 
shown in solid. At later times, the universe is homoge¬ 
neously reionized and the dashed line shows the kinetic 
SZ effect for this period. 

Inverse Compton scattering in galaxy clusters leads 
to a larger signal than that caused by Doppler scatter¬ 
ing. Because only a small fraction of the CMB photons 
present are likely to scatter inside the cluster medium 
(~ 1% for a massive 15 keV cluster), the distribution 
will not thermalize to that of the hot gas. The effect 
has a characteristic frequency dependence, and it can be 
distinguished from Doppler scatterings which leave the 
frequency distribution of the photons unchanged. The 
dot-dashed line in FigureElis the combination of the sec¬ 
ondary signal with the primordial CMB, once the ther¬ 
mal SZ with its characteristic frequency distribution has 
been removed. 

On scales where patchiness during reionization con¬ 
tributes more to the total signal than the Doppler effect 
owing to modulations in the density, {I < 4000), the pri¬ 
mordial anisotropies dominate. To study the secondary 
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Fig. 4.— Based on the implementation of Model B with a constant value of we divide the Doppler signal into three redshift epochs. 
Plot on the left shows the kinetic Sunyaev-Zel’dovich effect out to redshift 3, in the middle the Ostriker-Vishniac regime from z = 3 out to 
2 = 11 is shown, and on the right the regime of non-uniform (patchy) reionization. Each picture shows the same angular extent, 1°. The 
left panel exhibits highly clustered structures. With increasing redshift, the clustering becomes more linear and patchy reionization leads to 
CMB signals on large scales. We smoothed the maps obtained from our simulations by a Gaussian beam corresponding to I = 10000 ~ 1^ 
comparable to the angular resolutions of ACT and SPT. 



b/lO'*’ 


Fig. 5.— The probability distribution functions for the three 
maps in Figure]^ are shown. The solid line is for the kinetic SZ 
regime z=0-3, the long-dashed line for the Ostriker-Vishniac era 
z=3-ll, the short-dashed line for the patchy epoch of model B 
with constant ionization efficiency, z=ll-19. The two histograms 
representing the earlier cosmological epochs are well-approximated 
by Gaussians. 


anisotropies created during reionization one has to ob¬ 
serve angular scales where photon diffusion smoothes out 
the primary anisotropies. 

Upcoming experiments such as ACT and SPT will 
observe the sky in a number of frequency bands with 
comparable sens itivity and angular resolution. ACT 
llKosowskvll2003D will observe in three frequency bands: 
145, 225, and 265 GHz, with FWHM beam-widths 
&FWHM of 1.7, 1.1, and 0.9 arcminutes, respectively. We 
use the specifications of the 225 GHz channel for our 
power spectrum analysis because the thermal SZ effect 
almost vanishes in this frequency regime, its zero being 
at ~ 218 GHz. The sensitivity per resolution pixel 
for this channel is cr = 2 /iAT. The South Pole Telescope 
lIRiihl et a,lJl2'nnl will observe in five bands at 95, 150, 
219, 274, and 345 GHz. The 219 GHz channel will have 




Fig. 6.— Contribution to the kinetic SZ signal for Models A (up¬ 
per panel) and C (lower panel) out to different redshifts. Although 
a large fraction of the total Doppler effect comes from low redshift, 
high overdensity regions, there is a significant contribution out to 
2 ~ 10. In the patchy model (bottom) the regime z = 11 — 19 leads 
to further enhancement of the signal. 






















































Fig. 7.— Redshift dependence of the thermal Sunyaev-Zel’dovich 
AT/T = —2y effect as generated with our simulations. The ma¬ 
jor contribution comes from galaxy clusters in the regime 2 : < 3. 
Hence, the thermal SZ effect is less suited for studying the reion¬ 
ization epoch than is the kinetic SZ effect. 


an angular resolution of Ofwhm = 0.69' and a sensitiv¬ 
ity of lO/ritr ®. We assume sky coverages of 0.5% and 
10% for ACT and SPT, respectively. 

From these specifications and a template for the power 
spectrum of primary and secondary anisotropies at ar- 
cminute scales (which we obtain from our simulations), 
we can calculate the errors in the Ci determination, 
including noise as additional random held on the sky 
llTegmarldll997M : 


AC, = 


fsky{2l + 1) 


Cl 


fsky 


wBf 


(7) 


where we assume errors bars corresponding to a Gaus¬ 
sian map. Here, w = {0fwhmO')~^, and Bi is a beam 
prohle (assumed Gaussian) given by Bi = 

(with 01 , = 0f„jhm/V8 ln2). 

We assume that experiments such as ALMA ® will 
lead to a good understanding of the point source fre¬ 
quency spectrum and angular clustering. In more pes¬ 
simistic scenarios where only the frequency dependence 
or nothing is known about point sources, they can 
strongly degrade our ability to detect the kinetic SZ 
ijHiiffenberger fc Seliakl20d^ (we included their estimate 
for the IR source power spectrum at 145 GHz in Figure 
0. On the other hand the assumption of perfect cleaning 
of the thermal SZ effect is safer because we understand 
its frequency dependence well. 

To compare the experimental constraints with the 
power spectra extracted from our simulations on degree 
patches on the sky, we bin the errors into bands of width 
Al = 360. The predicted measurements errors are plot¬ 
ted in Figure 0 together with power spectra generated 
from the “cleaned” maps of primordial GMB and kinetic 


® J.E. Ruhl, private communication. 

® http://www.eso.org/projects/alma/science/alma-science.pdf 


SZ combined. Since the kinetic SZ appears almost fea¬ 
tureless on the scales accessible to the next generation of 
experiments, reionization models can be distinguished by 
their average amplitude on these scales. The error bars 
in the relevant band-powers of ACT are combined using 
^tot = (Ti I* follows that the overall amplitude 

can be measured with an accuracy of ctat = O.Oll/iAT 
(we assumed a 1% calibration uncertainty). Given that 
the plateau of the extended patchy model (Model C) lies 
at AT ~ 2.407fkK, while Model A has an average am¬ 
plitude on these scales of AT ~ 2.076/j,K, the two will 
be distinguishable at the 30 ct level with ACT, if we ig¬ 
nore contamination by point sources and by the thermal 
SZ, which is not a realistic assumption. This is the same 
for the South Pole Telescope, which has better angular 
resolution, but will look less deep (it has a much larger 
survey area). 

It is also of interest to ask whether additional con¬ 
straints on reionization scenarios could be achieved by 
using non-Gaussian statistics, in particular the four point 
function (the skewness arising from a line-of-sight veloc¬ 
ity effect should vanish). It may be expected, that be¬ 
cause the HII regions are created inside large overdensi¬ 
ties at high redshift, patchy models are more Gaussian 
on the whole. The main problem with this notion is that 
in any patchy model with a comoving bubble size compa¬ 
rable to ours this effect will suffer from “washing out” by 
the primordial GMB. On scales where Doppler induced 
fluctuations become larger than the primordial GMB, the 
impact of patchiness amounts only to a fraction of the 
total signal (15% in Model B, 30% in Model G), the rest 
being attributed to scatterings owing to density modu¬ 
lations alone. Also, the majority of non Gaussian con¬ 
tributions comes from clusters and filaments at 2 < 3, 
compare Figure 0 We computed the kurtosis 


04 


((AT/T)4) 

13 


( 8 ) 


of thermal SZ “cleaned” maps (i.e. primordial 
GMB-fDoppler) that were filtered with a Gaussian in 
Fourier space, approximately cutting out “contamina¬ 
tion” by the primordial GMB anisotropies below I ~ 
3500, and beam smearing at scales of I ~ 9000 compara¬ 
ble to the angular resoluti on of AGT (similar to the win - 
dow function proposed bv iHuffenberger fc Seliaklf2004ll . 
Using this method, we do not find a statistically signifi¬ 
cant difference between patchy and homogeneous models 
of reionization. 


5.2. Expected bias in cosmological parameter 
determination from Planck 

Gosmic microwave background measurements have 
been and will likely remain the most precise tools for 
the measurement of cosmological parameters. The best 
constraints will come from a combination of GMB tem¬ 
perature and polarization power spectra which encapsu¬ 
late all the relevant information in the sky maps. The 
aim is to make the data cosmic variance limited to as 
high as possible multipole numbers. WMAP is cosmic 
variance limited up to / ~ 500. 

The Planck satellite should achieve cosmic variance 
limitation out to I = 2500 for its temperature power 
spectrum measurement. Hence, Planck will reach into 
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Fig. 8.— Different contributions to the total Doppler signal (kinetic SZ) in our extended patchy reionization Model C. The solid curve is 
the contribution to temperature anisotropies from the patchy regime alone at z = 7 — 20. The dashed curve gives the Doppler effect from 
density modulations at homogeneous ionization out to z = 11. The dotted curve sums up those contributions to the total Doppler signal. 
The total temperature fluctuations, cleaned of the thermal SZ and IR sources (shown at 145 GHz in the thin dashed lines), is given by the 
dot-dashed curve. These curves are smoothed versions of the power spectra generated from 50 maps. 


the regime where the secondary anisotropies become im¬ 
portant. To avoid biases in parameter estimation, sys¬ 
tematic changes that the secondaries may produce in the 
power spectra need to be considered. 

The contribution to the power spectrum from patchy 
reionization is one to two orders of magnitude smaller 
than the primordial CMB anisotropies on the relevant 
scales, but owing to the exquisite sensitivity of these ex¬ 
periments it biases parameter estimates. For the analytic 
models of patchy reioniz ation suggested by iKnox et al.l 
(TT^ and lSantos et al ] this parameter bias was 

estimated. In these models the mean bubble size is 
smaller than in our computation, and the signal peaks 
at higher multipoles. 

The bias can be estimated from the Fisher matrix co¬ 
efficients 


I 


dCi dCi 

dpi dpj 


(9) 


in the following manner: 



CTi 




dCi 


cF 




( 10 ) 


where are the systematic biases in the determination 
of parameters pi, wi are the inverse squares of the sta¬ 
tistical errors in the Ci estimation, given by Equation 0 

and ai = \J~F^ is the estimate of the error bars for pa¬ 
rameter i. In this expression, CF denotes the combined 


model 

'^ri 

Ua 

^dm 

iOb 

ris 

As 

Icr 

0.0035 

0.010 

0.0017 

0.00018 

0.0045 

0.0050 

A 

0.40 

1.18 

-1.10 

1.71 

2.14 

0.16 

B 

1.26 

1.66 

-1.54 

2.40 

3.05 

0.78 

C 

2.25 

2.89 

-2.69 

4.20 

5.62 

1.58 


TABLE 1 

Bias in units of the statistical error [Bi) expected for 

COSMOLOGICAL PARAMETER ESTIMATION WITH PlANCK, IF 
temperature and POLARIZATION POWER SPECTRA ARE USED AND 
THE INFLUENCES OF kSZ/OV AND PATCHY REIONIZATION ARE 
NEGLECTED IN THE POWER SPECTRUM ANALYSIS. ThE MAXIMUM 
MULTIPOLE IN OUR ANALYSIS WAS I = 4000. 


Doppler power spectrum owing to kSZ/OV and patchy 
reionization. 

In Table ^ we show the results of our analysis, using 
the power spectra for our different models as bias 


The first line shows the statistical error alone, 



We combined the three frequency channels of Planck with 
the highest angular resolution and took into account the number 
of polarized instruments. For a one year observation period, the 
three channels (217, 143 and 100 GHz) have fffwhm = 5.0,7.! and 
9.2 arcminutes. This leads to the raw sensitivities: 

vj-^ = {0mS4)iKf 

Wp^ = {0.0200iiKf . 

We assumed a sky coverage of fsky = 0-8. 
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Fig. 9.— With future experiments like ACT, the various reionization models should be easily distinguishable by measuring the temperature 
power spectrum. Since the kinetic SZ signal is almost featureless on the scales of experimental relevance, the band-powers can be combined 
to a simple amplitude when distinguishing the different models. We find that the double reionization scenario (given by the upper power 
spectrum) could thus be distinguished at high significance from a uniform, instantaneous model where the universe reionizes at 2 = 14 
(bottom curve). The change of the ionizing efficiency with time C affects the slope of the secondary anisotropy power spectra. 


the following lines show the parameter bias Bi for Models 
A, B, and C. For Planck, we used the specifications of 
the High Frequency Instrument The power spectra 
derivatives for the Fisher analysis were computed for the 
fiducial cosmological model given in the introduction. 

It is clear from Table ^ that secondary anisotropies 
need to be taken into account. This is the case in par¬ 
ticular for parameters that influence the shape of the 
power spectra at intermediate scales (wf,, ^s)- Con¬ 

straints on the amplitude Ag and the optical depth owing 
to reionization Tri are less heavily biased (besides the am¬ 
plitude, the reionization optical depth only affects CMB 
polarization on large scales). Note that even in a sim¬ 
ple homogeneous reionization model, as our Model A, 
most of the biases are of order unity. The individual sys¬ 
tematic shifts depend strongly on what p arameters are 
used. The analysis of iSantos et alJ (1200,11) has extra pa¬ 
rameters that make the biasing source (the patchy power 
spectrum) distribute differently into each parameter off¬ 
set. 

The polarization anisotropies generated during the 
reionization epoch are expected to be four orders of mag¬ 
nitud e smaller than the temperature power spectra llHiJ 
l2?)00h . This suggests that to avoid biases one could use 
temperature information down to an angular scale where 
the Doppler contamination is still negligible but use the 
polarization data for all 1. We show the results of an 
analysis in which we used temperature information only 

http://www.rssd.esa.int/index.php?project=PLANCK&;page= 


model 

'^ri 

Ua 

^dm 

LJb 

ris 


Icr 

0.0041 

0.013 

0.0020 

0.00026 

0.0079 

0.0056 

A 

0.010 

-0.018 

0.017 

0.024 

0.086 

0.068 

B 

0.027 

-0.028 

0.026 

0.037 

0.135 

0.196 

C 

0.050 

-0.070 

0.067 

0.058 

0.301 

0.487 


TABLE 2 

Bias in cosmological parameter estimation with Planck is 
REDUCED significantly, WHEN TEMPERATURE POWER SPECTRA 
ARE USED ONLY UNTIL I = 1000 AND POLARIZATION POWER 
SPECTRA ARE USED IN THE WHOLE RANGE (PlANCK SHOULD BE 
ABLE TO MEASURE OUT TO I ~ 1800). On THE OTHER HAND, 

THE PARAMETER ERRORS BECOME ONLY SLIGHTLY LARGER BY 
LEAVING OUT THE HIGH I TEMPERATURE INFORMATION. 


out to I ~ 1000, but polarization in the full range acces¬ 
sible to Planck (this will be out to I ~ 1800) in Table 
12 The Icr error bars for the parameters are only slightly 
larger than in Tabled indicating that the parameter es¬ 
timates are more or less “saturated” at I ~ 1000. On 
the other hand, the expected parameter biases owing to 
reionization are much smaller. The use of polarization 
information in future CMB experiments thus can play 
an important role beyond breaking degeneracies between 
traditional cosmological parameters and improving the 
error bars. 

The other strategy is to include an extra parameter 
to model the effect of the Doppler contributions. This 
:perf.a(pproach has the added advantage that a positive con- 
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strain! on reionization could be obtained with Planck 
alone. Patchiness shifts the power spectrum on small 
scales, and we model this by a “patchy amplitude” pa¬ 
rameter. This parameter is included by adding the 
Doppler spectra given in the last section with a variable 
amplitude Ad io the primordial CMB. If the parameter 
derivatives in the Fisher matrix center {Ad = 1) around 
Model B, we find that the amplitude parameter could 
be constrained with gad = 0.51. The model is only 
~ 0.3(7 away from homogeneous reionization (Model A), 
so Planck will not be able to make a strong statement. 
If patchy reionization turns out to be extended, similar 
to Model C, Planck should observe this at higher sig¬ 
nificance a Ad = 0.20, amounting to a 5tT detection of 
Doppler induced secondary anisotropies, with the am¬ 
plitude of Model C being 3tT away from homogeneous 
reionization. Finally, if reionization proceeded homo¬ 
geneously, Planck will not gain knowledge beyond its 
large scale polarization measurements of by using 
small scale temperature fluctuations, given that the un¬ 
certainty for our Model A lies at aAo = 0.71. 

When Planck’s power spectra are modeled with the 
“patchy amplitude” parameter, constraints for the stan¬ 
dard cosmological parameters are slightly improved 
over the analysis that abandoned temperature data 
beyond I = 1000 (Table EJ, despite the introduc¬ 
tion of an extra parameter. Concretely, we find in 
this analysis that tT^^.=0.037, ctha =0.011, CT(^^^=0.0018, 
cr^^=0.00021, crn„, =0.0058, cta, =0.0052. 

The Doppler power spectrum can bias the result of 
parameter analyses with future CMB surveys. Planck 
may be able to improve constraints on models with an 
extended epoch of patchiness. We hope to have shown 
in this section that careful modeling of the epoch of first 
stars will become crucial for doing precision cosmology 
with progressively more refined CMB experiments. 

6. DISCUSSION 

We have presented simulations of secondary 
anisotropies in the cosmic microwave background 
calculated using smoothed particle hydrodynamics 
simulations of large scale structure, focusing on the 
effect produced by patchy reionization. We incorpo¬ 
rated an analytic model for the morphology of the HII 
regions into our numerical treatment and investigated 
whether such patchy scenarios can be distinguished from 
homogeneous reionization. 

An important advantage of our technique over pure an¬ 
alytical predictions of patchy reionization morphology is 
that we follow the complicated clustering of dark matter 
and baryons into the slightly non-linear regime. In con¬ 
trast to full radiative transfer calculations of the reion¬ 
ization epoch, we combine uncertainties in the physics 
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of source formation, feedback processes and radiative 
transfer into a single parameter and explore the conse¬ 
quences of varying this parameter. This simplification 
allows us to make predictions on scales an order of mag¬ 
nitude larger than current radiative transfer schemes can 
accomplish with a small expense of memory and CPU. 

We extracted power spectra from sky maps produced 
by tracing rays across our simulation volumes. The 
patchy reionization signal peaks on multipole scales of 
I ~ 2000, and it increases the amplitude of the “cleaned” 
CMB power spectrum by up to 30% on scales I > 4000, 
so that the total level of Doppler related anisotropies 
is AT ~ 2.4/iAr. We found that with the next gener¬ 
ation of ground based CMB experiments (ACT, SPT) 
the different reionization models we investigated could be 
distinguished with high significance by using the power 
spectrum. Additional information about the morpholog¬ 
ical properties of that epoch that could be obtained for 
instance by measuring the four point function or other 
deviations from Gaussianity will probably be difficult to 
obtain, because the reionization signal peaks on angular 
scales where the primordial CMB anisotropies dominate. 
In the future, it may be possible to combine measure¬ 
ments of the CMB with other observat ions such as 21 cm 
fluctuations from neutral hydrogen llZaldarriaga et ^ 
l2?inl ICoora,l[2f)nl o r Lyman-alpha emission from high 
redshift galaxies (e.g. iFurlanetto et al.lf2004(Jl to further 
constrain the topology of the reionization process. 

We investigated the bias in the determination of cos¬ 
mological parameters that will be produced by the addi¬ 
tional patchy reionization signal, when extracting cosmo¬ 
logical parameters from CMB anisotropy measurements. 
We find that for Planck this bias is significant. The bias 
may be circumvented by focusing completely on polariza¬ 
tion information in the multipole regime where patchi¬ 
ness peaks, with only a slight disadvantage in parameter 
constraints. Alternatively, a template for the Doppler 
spectrum could be introduced in the parameter analy¬ 
sis which may lead to a detection of the effects of an 
extended reionization phase. 
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